Tuning magnetic avalanches in Mni2-ac 
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Using micron-sized Hall sensor arrays to obtain time-resolved measurements of the local magneti- 
zation, we report a systematic study in the molecular magnet Mni2-acetate of magnetic avalanches 
controUably triggered in different fixed external magnetic fields and for different values of the initial 
magnetization. The speeds of propagation of the spin-reversal fronts are in good overall agreement 
with the theory of magnetic deflagration of Garanin and Chudnovsky [Ij]. 
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1. INTRODUCTION 

Mni2-acetate is a crystal composed of magnetic 
molecules, each of which behaves as a high-spin, high- 
anisotropy magnet [1, Q. At low temperatures, the 12 
Mn atoms are strongly coupled via superexchange to 
form a ferrimagnet with a net (rigid) spin 5 = 10. Mag- 
netic interactions between the molecules are thought to 
be negligible so that Mni2 can be modeled by an effective 
spin Hamiltonian: 



n = -DSl - ASt - gtMBS^B, 



Ha 



(1) 



where Bz is a magnetic field applied along the c-axis of 
the crystal, Sz is the z component of the spin, D — 0.548 
K, A= 1.17 X 10~^ K, g = 1.94, and Ti± represents small 
symmetry-breaking terms that allow tunneling across the 
anisotropy barrier 0, [E, @] . The energy barrier against 
magnetic reversal, U, is easily calculated from Eq. [T] 
The magnetic relaxation rate for an individual molecule 
becomes sufficiently slow at low temperatures (< 2 K) 
that the magnetization of the crystal can be prepared and 
maintained in a metastable state for time periods well in 
excess of the experimental times. Once in this metastable 
state, the magnetization may relax as an abrupt (< 1ms) 
"magnetic avalanche," an exothermic process involving 
the release of Zeeman energy Q . These spatially inhomo- 
geneous reversals proceed as a traveling "front" between 
regions in the crystal with opposing magnetization and 
have been described as magnetic deflagration in analogy 
with chemical deflagration [l|, Q . 



In general, the speed of a deflagration front is gov- 
erned by two parameters: the thermal diffusivity, k, 
which specifies the rate at which heat diffuses through 
the medium, and the reaction rate of the constituents, 
r(T/), where T/ is the "flame temperature" produced 
by the reactants near the front. Combining these pa- 
rameters gives an approximate expression for the speed, 
V ~ -s/kT 0. In the case of magnetic deflagration, the 
medium through which heat flows is the crystal, and T is 
the relaxation rate of the metastable spins, which obeys 
an Arrhenius law. 



r = ro exp [-U{B)/T].^ 



(2) 
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with To = 3.6 X lO'^ s'^ The relaxation rate can 

be increased by lowering the barrier with an increasing 
external magnetic field or by increasing the temperature, 
T. Although K has not been measured for Mni2 -ac, a 
value of K ~ 10^^ m^/s was deduced from the avalanche 
data in Ref. [l3| ■ Suzuki et al. [8| showed that the speeds 
of magnetic avalanches can be modeled approximately as, 
V ^ a/kTo exp[— [//(2T/)], where T/ is the temperature 
at or near the propagating front associated with the en- 
ergy released by the reversing spins. The theory of mag- 
netic deflagration stands in qualitative agreement with 
experiments, yet more precise quantitative confirmation 
remains an open experimental challenge [l|, d, [l^, [13, , 
which is undertaken here. 

There are two parameters under experimental control 
when performing magnetic avalanche studies on Mni2-ac: 
the external magnetic field and the initial magnetization 
which tunes the metastable spin density and thus the 
available "fuel" . Varying the external magnetic field af- 
fects both the barrier against spin reversal and the energy 
released which determines the temperature Tf produced 
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by the reversing spins. The variation of the metastable 
spin density affects primarily the energy released and 
therefore, Tf. By varying these parameters indepen- 
dently, we are able to explore the wide ranging conditions 
in which avalanches may be triggered. In particular, we 
report systematic studies for three classes of avalanche 
preparations. The first class (I) contains avalanches trig- 
gered at various external fields with fixed (maximum) 
initial magnetization. For these avalanches, both U and 
Tf vary. The second class (II) contains avalanches trig- 
gered at a fixed external field, but with various initial 
magnetizations. Avalanches of this class differ primarily 
through Tf, with U varying only through the internal 
fields. And the final class (III) contains avalanches trig- 
gered with various initial magnetizations and external 
fields, such that Tf is approximately constant while U 
varies (details below). 

We report the speeds of propagation of the fronts for 
these three classes of avalanches, allowing us to make a 
thorough comparison with the theory of magnetic defia- 
gration. Our results are in overall agreement with the 
theory. With certain simplifying assumptions detailed 
below, we obtain temperatures between 5 and 16 K and 
thermal diffusivities ranging from k — 1.2 x 10~^ to 
8.7 X 10~^ m^/s, consistent with the value of k estimated 
inRef. [3. 
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FIG. 1: The initial magnetization, M^, as a function of ap- 
plied cooling field. Msat is defined as the magnetization with 
all the spins aligned in the positive direction. The change 
of magnetization during an avalanche is AAf = Alsat — M^. 
The same curve is found for all crystals measured. The inset 
shows the energy levels biased with a -1-0.3 T field. 



2. EXPERIMENT 

The magnetization dynamics were studied with an ar- 
ray of 30 X 30 /im^ Hall sensors spaced apart by 80 /im 



center-to-center [T5|. Crystals of Mni2-acetate with di- 
mensions about 1.0 X 0.2 x 0.2 mm"^ were attached to the 
array with Apiezon M grease. The crystal was encased in 
grease along with a constantan wire placed near the sam- 
ple for use as a heater in an arrangement similar to that 
used in Ref. [T^ . The entire assembly was immersed in 
•^He at temperatures down to < 300 mK. 

Prior to triggering an avalanche, the sample was pre- 
pared in a metastable magnetic state. To do this, the 
sample was cooled from a high temperature (~ 6 K) down 
to 300 mK in the presence of a small external magnetic 
"coofing" field between ±0.3 T. The inset of Fig. [1] is 
a schematic of the energy levels biased with a -1-0.3 T 
cooling field; as shown in Fig. [U the magnetization of 
the sample, M^, depends on the magnetic field in which 
it was cooled. At 0.3 K, only the Sz = ±10 states are 
appreciably occupied. Mz reflects the ratio of these oc- 
cupied states. Once the sample is well below the blocking 
temperature, the external field can be changed without 
changing the magnetization (the system is blocked). The 
field is then increased to a predetermined value (> +1.25 
T) [g. When the field has stabilized (~ 1 min), a cur- 
rent is passed through the wire heater gradually raising 
the temperature and triggering the avalanche. For more 
details on triggering avalanches with this method, see 
Ref. [14]. 

All avalanches reported here were triggered in a posi- 
tive field. The amount of metastable magnetization that 
reverses during an avalanche is given by AM = \Msat — 
Mz\, where Mz is the initial magnetization and Mgat is 
the magnetization with all spins aligned in the positive 
direction. For full magnetization reversal, AM = 2Msat- 
For convenience, we introduce the parameter AM /2Msat 
as a dimensionless measure of the initial metastable mag- 
netization density. As an example, cooling the sample in 
zero field leads to Mz = 0, or AMj^Msat = 0.5. 

We prepared the three classes of avalanches by varying 
the initial magnetization, M^, and the external magnetic 
field Hz- Variation of Hz tunes the barrier, [/, as well 
as the average energy released per molecule during an 
avalanche, 

where Bz = lJ.o{Hz + Mz). Avalanches of class I are 
those with initial magnetization Mz = —Msat, i-e., 
AM/2Msat — 1, and Hz is varied. Class II avalanches 
are triggered at a fixed external field, with various ini- 
tial magnetizations. Finally, class III avalanches are trig- 
gered at various Hz and Mz, such that (E) remains con- 
stant. 

We collected avalanche data on four different crystals 
with dimensions 1.00 x 0.20 x 0.20 mni"^ (Sample A), 
1.20x0. 10x0. lOmm^ (Sample B), 0.80 x 0.15 x 0.15 mm^ 
(Sample C), and 1.00 x 0.25 x 0.25 mm^ (Sample D). We 
report detailed data on crystal A. Although the absolute 
values of the avalanche speeds differed for different crys- 
tals (as discussed in detail later in this paper), similar 
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FIG. 2: The time response of five equally spaced Hall sensors 
placed along the length of the crystal. The peak indicates 
the arrival of the magnetization interface (avalanche front). 
The inset shows the traveling front that separates spin-up 
and spin-down regions. Note the transverse magnetic field 
lines at the position of the spin-reversal front. 



behavior was obtained for all crystals as a function of 
the experimental parameters. 



3. RESULTS 

As was shown by Suzuki et al. [1] , the avalanche pro- 
gresses through the crystal in a similar fashion to that of 
a domain wall in a ferromagnet (l7j . There is an inter- 
face separating regions of opposing magnetization, which 
produces a large transverse magnetic field, near the 
front, as shown schematically in the inset of Fig. O Fig- 
ure [2] shows the time response of five equally spaced Hall 
sensors produced by a zero field cooled avalanche. At 
time t = 0, the magnetization is zero. At t « 0.32 ms, 
the magnetization begins to reverse near the first sen- 
sor as the avalanche front approaches. At t = 0.41 ms, 
the signal on the first sensor is maximum indicating the 
avalanche arrival at position fun. By t — 0.65 ms all 
spins have reversed and the sample is completely magne- 
tized. The avalanche speed is deduced from the arrival 
time of the peak at each sensor and the known spacing 
between the sensors. 

Figure [3] shows the avalanche speed for class I 
avalanches triggered at various fields, all with the same 
initial metastable magnetization, AM/2Msat — 1- The 
speeds of the avalanches increase as the field is increased. 
This is expected since increasing the field both lowers the 
barrier and increases (E). There are local maxima at 
fields corresponding to the tunneling resonances, which 



FIG. 3: Propagation speed as a function of applied magnetic 
field for Class I avalanches in crystal A for Class I avalanches, 
AM/2Msat = 1 
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FIG. 4: Propagation speed as a function of AAI /2Msat for 
Class II avalanches in crystal A. The field was fixed at 2.50 T 
and 2.20 T, while AM/2Msat was varied between about 0.10 
and 1.00. 



are denoted with vertical dotted lines. This tunneling 
enhancement of the avalanche speed is consistent with 
previously reported results p,2, 13,, JJJ. It is also an in- 
dication that the flame temperatures are low enough to 
preserve the effective Hamiltonian, Eq. [1] and the rigid 
spin (S = 10) approximation. 

Figure [5] shows the speed for class II avalanches trig- 
gered at fixed values of the external field. In particular, 
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FIG. 5: Propagation speed as a function of applied magnetic 
field for Class III avalanches in crystal A. The initial magne- 
tization and external field are adjusted to hold the average 
energy released per molecule constant at (iJ) = 32.6 K and 
18.3 K. The vertical lines indicate the fields at which quantum 
tunneling occurs for Mni2-ac. Note that the avalanche speed 
displays clear oscillations as a function of magnetic field, with 
higher values on-resonance than off-resonance due to quantum 
tunneling. This can also be seen in Fig. [T] 
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FIG. 6: Speed of propagation of Class II avalanches triggered 
in four different crystals at fioH^ = 2.5 T for various initial 
magnetizations. 



FIG. 7: Avalanche speed vs. the ratio U/Tmax, where U is the 
calculated barrier and Tmax is the maximum fiame tempera- 
ture. The inset shows Tmax calculated from the heat capacity 
at 2.5 T as a function of (E). 



data are shown for fioH^ = 2.50 T and iioHz = 2.20 T. 
By keeping the field fixed and varying the initial magne- 
tization, only (E) varies while U remains approximately 
fixed. As AM/2Msat decreases, so does the speed. 

Figure [5] shows the speed for class III avalanches trig- 
gered at various external fields. The initial magnetiza- 
tions were also varied such that (E) remained approx- 
imately constant. Presumably, the flame temperature 
is also nearly constant for all avalanches so prepared. 
Therefore, the variation in avalanche speeds should be 
due to variation in the field-dependent barrier U . Again, 
the vertical dotted lines drawn on Fig. [5] denote the lo- 
cation of the resonant fields for Mni2-ac. 

Figure [6] compares the class II avalanche speeds for all 
four crystals. The speed of the avalanche varies consid- 
erably from sample to sample; however, any dependence 
on the sample dimensions is not immediately obvious. 



4. COMPARISON WITH THEORY AND 
DISCUSSION. 

Garanin and Chudnovsky [l| developed a comprehen- 
sive theory of magnetic deflagration describing the ig- 
nition and propagation of the deflagration front. For a 
planar front, the avalanche speed is given over a broad 
range of parameters by the simple approximate expres- 
sion [11: 



l3KTfr{B,Tf) 



U{B) 



(4) 
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FIG. 8: (a) Avalanche speed vs. the ratio {U* /Tj), where the scaled barrier U* — {1 — oi-j^)U with a = 0.13, and the 
flame temperature Tf — 0.67 x r,„ai with k — 1.2 x 10^^ m^/s. The inset illustrates that a data collapse can be obtained 
by adjusting [/*, (here k = 2.4 x 10~^ m^/s); however, numerical agreement with the measured data requires the additional 
scaling parameter Tf. (b) Data collapse obtained using a calculated from the measured transverse field during the avalanche, 
setting T = Tmax, and allowing the thermal diffusivity to vary with temperature. The fit yields k = (1.8 x 10^*) T^'^ m'^/s. 



where Tf is the temperature at the front (flame tem- 
perature), K is the thermal diffusivity, U{B) is the field- 
dependent barrier in units of Kelvin, and T{B, Tf) is the 
relaxation rate of the mctastable spins (Eq. [5]). 

Ref. [l^ established the importance of the dipolar 
fields in Mni fully magnetized sample adds (or 

subtracts) ^qM^ = ±52 mT to the applied external field, 
fioHz- We account for the initial magnetization, M^, 
when calculating the various field dependent quantities 
using Bz = noiHz + IVh), where -52 mT < < 52 
mT (see Fig. [l]). The barrier U{Bz = fioH^ + ^iqM^) is 
calculated from the effective spin Hamiltonian (Eq. [1]). 

The average magnetic energy released by the relaxing 
spins (Eq. [3]) leads to an increase in the temperature near 
the front. Assuming no heat loss, the maximum possible 
temperature, Tmax, can be calculated using the experi- 
mental heat capacity reported in Ref. [ll| . The heat ca- 
pacity of Mni2-ac depends on the magnetic field. There- 
fore, we subtract the calculated zero-field spin (Schottky) 
contribution from the measured zero-field heat capacity 
from the data reported in Ref. [HI. To this we add the 
calculated spin contribution at a specified field, B^, for 
the total field dependent heat capacity Ctot{Bz,T). By 
equating the integral of this heat capacity to the aver- 
age energy released per molecule, (i?), we can calculate 



{E) = / Ctot{Bz,T)dT. 
Jo 



(5) 



We assume the initial (ignition) temperature is much less 



than Tmax- This is a reasonable approximation, as the 
ignition temperatures for avalanches triggered above IT 
are below 1 K (see ref. [l3|) compared with values calcu- 
lated for Tmax between 7 and 18 K (depending on (E)). 
Tmax (for fioHz = 2.5 T) is shown as a function of (E) 
in the inset of Fig. [7] (a). 

We now proceed to compare our data with the theory 
of Garanin and Chudnovsky py] as given by Eq. 21 If we 
assume that the thermal diffusivity k is a constant (or 
a weak function of temperature), then the speeds for all 
avalanches should lie on a single curve when plotted as 
a function of U/Tmax- Figure [7] shows avalanche speeds 
for crystal A for the three different experimental proto- 
cols shown in Figs. [31 [H and O plotted as a function 
of U/Tmax- Although the overall behavior for the three 
types of avalanches is similar, the data do not lie on one 
curve. 

The deviations could arise from several factors. (1) We 
have used an Arrhenius form for the magnetic relaxation; 
departures from Arrhenius law behavior are unlikely to 
be responsible for the deviations as it has been found 
experimentally to hold reasonably well in the range of 
temperature of our experiment [lol . [ll| . (2) The thermal 
diffusivity is known to depend on temperature, while we 
have assumed it to be constant. (3) The fiame temper- 
ature may be lower than the value calculated from the 
specific heat, as some of the energy may escape the sam- 
ple, or be distributed ahead of the front. (4) The barrier 
U may be reduced by the transverse component of the in- 
homogeneous field, Bx, established at the traveling front 
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by the reversing spins (see inset to Fig. [2]). The effects 
of can be included in the calculation of U by includ- 
ing an additional Zeeman term {—g^iBSxBx) in Eq. [1] 
In addition, provides a symmetry-breaking term that 
increases the tunneling rate [13, [2l| . 

We note that the deviations shown in Fig. [7| are espe- 
cially pronounced at high values of U/Tmax- In particu- 
lar, the class I avalanches with AM/2Msat = 1 (shown 
as filled diamonds) have the highest speeds. This sug- 
gests that reduction of the potential barrier U for large 
AM/2Msat plays an important role. 

It is unclear how to incorporate the effects of a spatially 
inhomogeneous transverse field component into the ana- 
lytical theory of magnetic deflagration (Eq. |4]). Instead, 
we include the effects of Bx on the relaxation rate by 
introducing an effective barrier: 



U* 



AM 



f/, 



(6) 



where a is determined empirically ^22j . Although the 
scaling factor aAM/2Msat explicitly contains Mz (thus 
appearing to be used twice), it is used here only to ac- 
count for the size of Bx- This is justified by our exper- 
iment, since the maximum value of measured by the 
Hall sensors during an avalanche is found to be propor- 
tional to AM/2Msat- 

The inset to Fig. [H] (a) demonstrates that a collapse 
onto a single curve is obtained for a = 0.13 ±0.01. How- 
ever, the collapsed curve does not agree with the theory, 
shown by the dashed curve. An additional step can bring 
them into line, as described below. 

As pointed out earlier, due to possible heat loss 
through the edges of the crystal and/or heat diffusion 
ahead of the front, the flame temperature Tf may well 
be less than T^ax- Assuming that Tf is proportional to 
Tmax, the constant of proportionality is deduced from fit- 
ting the data in the inset of Fig. [8] (a) with Eq. [H The 
diffusivity, still assumed to be temperature- independent, 
is also treated as a fitting parameter. As shown in the 
main part of Fig. [8] (a), agreement with theory is ob- 
tained for crystal A for Tf = (0.67 ± 0.02) x Tmax and 
K = 1.2 X 10~^ m^/s . Using this analysis, all crystals 
show similar dependence on the barrier U* and the flame 
temperature Tf and yield Tf « 0.67 x Tmax- However, we 
find that the thermal diffusivity ranges from 1.2 xlO~^ 
to 8.7 xlO~^ m^/s from crystal to crystal. 

The fact that Tf is the same fraction, 0.67 x Tmax, 
for all crystals is a puzzle. The rate at which heat es- 
capes the crystal during an avalanche must affect Tf- 
The rate of heat loss is controlled primarily by the crystal 
cross section, surface roughness, and the thermal mount- 
ing conditions. Variations in the mounting conditions 
inevitably occur (e.g., thickness of insulating grease), al- 
though every effort was made to use similar conditions 
from crystal to crystal. There were no obvious visible dif- 
ferences in surface quality of the crystals. The cross sec- 
tions, however, were deliberately varied from 0.10 x 0.10 
to 0.25 X 0.25 mm^. One expects that the crystals with 



smaller cross sections should lose more heat through the 
boundaries and should consequently have lower flame 
temperatures and, according to Eq. [H smaller speeds. 
Figure [5] shows that the maximum speeds vary by ap- 
proximately a factor of 2.5 from crystal to crystal, but 
without the expected dependence on cross section. This 
implies that the widely different avalanche speeds in the 
four crystals (see Fig. 6) are unlikely to be due primarily 
to heat loss. Instead, we suggest that the variation of the 
avalanche speeds are attributable to variations of k. The 
thermal diffusivity of dielectric crystals (like Mni2-ac) at 
low temperatures is known to be strongly dependent on 
the defects, surface roughness, and dislocations in the 
crystal 

An additional puzzle is the large amount by which the 
potential barrier needs to be reduced to obtain a fit by 
the above analysis. From our data, we deduce a bar- 
rier U* that is 87% of the classically calculated barrier, 
U- A straightforward calculation implies that a trans- 
verse field of « 0.4 T is required to reduce the barrier by 
that amount. The largest Bx field recorded by the Hall 
sensors during an avalanche is only ~ 0.05 T, an order 
of magnitude smaller. Although it may contribute to it, 
the measured transverse field cannot by itself account for 
the large reduction of the barrier. 

In the analysis presented above, the thermal diffusiv- 
ity was assumed to be independent of temperature. We 
now relax this condition. We assume that the barrier U* 
is reduced by a much smaller amount corresponding to 
the measured transverse field, we set Tf = Tmax, and 
we allow K to assume a temperature dependence that 
yields the best fit. The result of this alternate fitting 
procedure, shown in Fig. [5](b), yields a collapse that is 
acceptable within the experimental uncertainties of the 
measurements. 

Remarkably, the latter method of analysis yields a 
thermal diffusivity that increases with increasing tem- 
perature approximately as k cx T'^. This form seems 
quite unphysical, as the thermal diffusivity normally de- 
creases as the temperature is raised. We suggest that 
this unexpected behavior may be associated with a spin- 
phonon bottleneck. 

A number of experiments have provided evidence that 
a spin-phonon bottleneck strongly affects the spin dy- 
namics and energy relaxation at low temperatures in 
molecular magnets such as V15 24 1, Feg [H, |2^ and Ni4 
[13] ■ In this process, the Zeeman energy generated by 
the reversing spins does not find a sufficient number of 
phonon modes at low temperatures to allow direct energy 
relaxation and equilibration, so that thermal equilibrium 
is established slowly while the energy is "bottlenecked" in 
the spin system. This bottleneck is lifted as the temper- 
ature increases, so that the number of available phonon 
modes increases and the energy is able to relax by direct 
spin-phonon processes. The effect of the bottleneck can 
find expression within our analysis as either a departure 
from Arrhenius Law behavior (which we have assumed to 
be valid) , or as an anomalous temperature dependence of 
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the thermal difFusivity. Within this scenario, k would not 
display the same temperature dependence when obtained 
by the standard method of measurement where one de- 
termines the time of propagation of a heat pulse, since in 
this case the energy is deposited into the phonon system 
directly. 

5. CONCLUSION 

We have presented the results of a thorough investiga- 
tion in Mni2-ac of the behavior of magnetic avalanches - 
the rapid reversal of magnetization that spreads through 
the crystal at subsonic speeds as a narrow interface be- 
tween regions of opposite spin. A controlled set of mea- 
surements in which some parameters were held fixed 
while others were varied provided systematic informa- 
tion, enabling a rigorous comparison with the theory. 

Two different methods were applied to fit data to the 
theory of magnetic deflagration of Garanin and Chud- 
novsky In the first, we suggest that the internal 
transverse field produced by the avalanche front affects 
the speed of the front itself. We model this effect with 
a reduced barrier, J7*, that varies with the size of the 
transverse field, Bx, produced by the avalanche front. 
Assuming a constant, temperature- independent thermal 
diffusivity, a reduced barrier U* allows a collapse of all 
the data onto a single curve. However, the transverse 
field measured at the front is not sufficiently strong to 
account for the large barrier reduction needed to obtain 
a good fit. A barrier reduction of this magnitude, if cor- 
rect, defies a simple classical analysis and may be a signal 
that quantum effects are important in the deflagration 
process for all values of fioH^, not just those associated 



with the tunneling resonances [1, [T^, [H, [31 ■ 

An alternative method of analysis that assumes a 
smaller barrier reduction commensurate with the mea- 
sured values of transverse field yields a temperature de- 
pendent K (X T^. We speculate that this rather surpris- 
ing temperature dependence may be real and due to a 
phonon bottleneck that becomes less effective as the tem- 
perature is raised. 

To summarize, we find overall agreement between our 
measurements and the theory of magnetic deflagration of 
Garanin and Chudnovsky However, detailed compar- 
ison yields either (A) a stronger reduction of the potential 
barrier than can be justified by the measured transverse 
fields; or (B) a thermal diffusivity that unexpectedly in- 
creases with increasing temperature, perhaps due to a 
phonon bottleneck; or (C) a combination of these (and 
possibly other) factors. Further confirmation of the the- 
ory and a better understanding of the avalanche process, 
would be provided by a detailed theoretical analysis of 
bottleneck effects, and independent measurements of var- 
ious parameters such as the flame temperature and the 
thermal diffusivity. 
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